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Abstract. We calculate the formation rates of low-mass 
X-ray binaries with a black hole. Both a semi-analytic and 
a more detailed model predict formation rates two orders 
of magnitude lower than derived from the observations. 
Solution of this conundrum requires either that stars with 
masses less than 20 Mq can evolve into a black hole, or 
that stellar wind from a member of a binary is accom¬ 
panied by a much larger loss of angular momentum than 
hitherto assumed. 
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1. Introduction 

Six low-mass X-ray binaries currently are known to have a 
mass function that indicates a mass of the compact object 
larger than ~ SMq, and these are therefore believed to 
have a black hole accretor. All these systems are soft X- 
ray transients and the total number in the galaxy of such 
systems is estimated to be between a few hundred and a 
few thousand; at any moment in time only a small fraction 
is X-ray active (for reviews, see Cowley 1992, Tanaka & 
Lewin 1995, Tanaka & Shibazaki 1996). 

The most popular formation scenario for low-mass X- 
ray binaries proposes a relatively wide binary with an ex¬ 
treme mass ratio as the progenitor system (van den Heuvel 
1983). The massive star evolves to fill its Roche lobe and 
engulfs its low-mass companion. Due to the extreme mass 
ratio the low-mass star spirals into the envelope of the 
high-mass star. A close binary remains if the spiral-in 
ceases before the low-mass companion coalesces with the 
compact helium core of the primary. The helium core con¬ 
tinues its evolution and may turn into a neutron star or 
a black hole. Only a few studies apply this evolutionary 
scenario specifically to the formation of low-mass X-ray 
binaries with an accreting black hole (see Romani 1992). 
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In this paper we discuss some problems of this stan¬ 
dard formation scenario in its simplest form (Sect. 2) 
and show that these persist in a more refined treatment 
(Sect. 3). We quantify this in Sect. 4 and in Sect. 5 the 
results are discussed. 

2. Semi-analytic approach 

In the standard scenario the formation of a low-mass X- 
ray binary starts with a high-mass primary and a low mass 
secondary star {rrio ^ ^Mq) in a detached binary. For 
simplicity we assume that the initial orbit is circular. As 
the primary evolves and expands to giant dimensions it 
fills its Roche lobe and starts to transfer mass onto its 
companion. Due to the small mass ratio mass transfer is 
highly unstable and results in a common-envelope phase. 
A close binary remains provided that the primary’s en¬ 
velope is fully ejected before the secondary star coalesces 
with the compact core of the primary. The reduction in the 
semi-major axis during the spiral-in can be computed by 
comparing the binding energy of the primary’s envelope 
with the orbital binding-energy of the binary (Webbink 
1984). The helium core continues its evolution and finally 
collapses into a neutron star or a black hole. The binary 
becomes an X-ray binary once the secondary star hlls its 
Roche lobe. 

We illustrate this scenario for primary stars of initially 
Mo = 20 and 60 M0which reach a maximum radius of 
1000 i?0(see Romani 1992). The scenario requires that the 
primary fills its Roche lobe during its evolution, i.e. that 
its Roche lobe has a radius less than 1000 Rq. The cor¬ 
responding semi-major axis of the binary orbit may be 
calculated using the equation for the radius Rl of the 
Roche lobe given by Eggleton (1983): 

Rl _ j, , ^ _ 0.49 

a " '^^‘^^"0.6 + g2/3in(l + g~i/3)- 

Here a is the semi-major axis of the binary and q = m/M 
is the mass ratio and Rl the Roche-lobe radius for the 
star with mass M. For the 20 and 60 M0primaries we 
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thus find a semi-major axis of about 1590 and 1440 Rq, 
respectively. 

The scenario further requires that the binary sur- 
vives the spiral-in that follows upon the first Roche- ^ ^ / 

lobe contact, i.e. that both the helium core and the ^ 

1 M 0 companion star are smaller than their Roche lobes. 

The mass and radius i?c of the helium core of the pri- 5 
mary can be computed with (see Iben & Tutukov 1985 
and de Loore & Doom 1992, respectively) 


Me = O.OTSM^-^^ (2) 

and 

logi?c = -1.13-b 2.26 log Me-0.78(logMc)2. (3) 

The mass and radius of the secondary star are not affected 
by the spiral-in. The mass ratio after the spiral-in gives 
the sizes of the Roche lobes of the helium core and the 
secondary in units of the semi-major axis. From the radii 
of the helium star and its companion we can thus derive 
the minimum semi-major axis for which both stars fit in¬ 
side their Roche lobes. For the 20 and 60 MQprimaries 
we thus find a semi-major axis after the spiral-in which 
should exceed 4.0 and 6.3 i?©, respectively. In both cases 
the most stringent limit is set by the main-sequence star. 

The ratio of the semi-major axes before and after the 
spiral-in can be computed by comparing the binding en¬ 
ergy of the primary’s envelope with the binding energy 
released by the shrinking binary (see Webbink 1984): 

f 1+ Me \ 

ai M \ aXR m J 

Here at and a/ are the semi-major axes at the onset and 
end of the spiral-in, and M and R are the mass and radius 
of the primary at the onset of spiral-in. We set aX = 0.5. 
We assume that the primary did not lose any mass before 
it fills its Roche lobe: M = Mo and Me = Mq — Me. 
For the 20 and 60 MQprimaries we thus find a semi-major 
axis before the spiral-in which should exceed 1470 and 
3160 Rq, respectively. 

For the 20 MQprimary, the lower limit on the semi¬ 
major axis for survival of the spiral-in is smaller than the 
upper limit for Roche-lobe contact and a low-mass X-ray 
binary can be formed when the initial semi-major axis is 
between these two limits. For the 60 MQprimary, the lower 
limit on the semi-major axis for survival of the spiral-in 
is larger than the upper limit for Roche-lobe contact and 
no initial orbit leads to the formation of a low-mass X-ray 
binary: Roche-lobe overflow always leads to a merger. 

Figure 1 shows, as function of the mass of the primary, 
the lower and upper limits to the semi-major axis of the 
initial binary at which the binary survives the spiral-in 
and the primary reaches the Roche lobe. Only for low 
mass primaries (M ^ 22Mq) does the binary survive the 
spiral-in. Thus, if only stars with an initial mass larger 
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Fig. 1. Lower limit to the initial semi-major axis at which the 
binary survives the spiral-in, as function of the mass of the 
primary. The upper (lower) dashed line gives the limit deter¬ 
mined from the condition that the secondary star (helium core) 
is smaller than its Roche-lobe, The solid line gives the upper 
limit at which the primary reaches its Roche lobe at its maxi¬ 
mum radius of Rmax = lOOOR©. The secondary is assumed to 
be a 1 M 0 star 

than 40 M©(van den Heuvel & Habets 1984) form a black 
hole, the formation of a low-mass X-ray binary with a 
black hole as a compact object is excluded in the standard 
scenario. 

One may argue that the binary survives the spiral-in 
even if the main-sequence star is larger than its Roche 
lobe after spiral-in, as long as the helium star fits inside 
its Roche lobe . The main-sequence star can shrink within 
its Roche lobe by transferring mass to the helium core. 
This mass transfer is stable because the main-sequence 
star is less massive than the helium core. In Fig. 1 we 
also show the lower limit to the semi-major axis obtained 
from the condition that only the helium star fits inside 
its Roche lobe. We see that even in this case low-mass 
X-ray binaries with a black hole cannot be formed in the 
standard scenario. 

3. Detailed model 

Since the stellar parameters used in the previous Sect, are 
rather rough, the same computation is performed with a 
more detailed model for population I (Z = 0.02) stars. 
We use the models with moderate core overshooting com¬ 
puted by Schaller et al. (1992), using the radius (calcu¬ 
lated from effective temperature and luminosity) and the 
mass of the star in the tabulated points. The mass of the 
star decreases as a function of time due to stellar wind. 
The mass loss in the stellar wind causes an increase of the 
Roche-lobe radius of the primary, (mainly) by increasing 
the semi-major axis and (to a lesser extent) by increasing 
the mass-ratio m/M. The increase in the semi-major axis 
is described, assuming an isotropic wind with high veloc- 
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itv according to the Jeans approximation, with (van den 
Heuvel 1983): 

a{M + m) = constant. (5) 


The relation between the semi-major axis at which the 
primary fills its Roche lobe Rl and the initial semi-major 
axis Oo is thus given by: 


Rl = RL{q)ai 


Rl 


m N Mo + m 
m) M + m' 


( 6 ) 


Here Mo and M are the zero-age mass of the primary 
and its mass at the moment it fills its Roche lobe. For 
each tabulated point of the evolutionary track we equate 
the radius R of the primary to the Roche-lobe radius Rl 
in Eq. 6 to calculate the corresponding maximum initial 
semi-major axis Oq. The values for Uo and R for a 20 and a 
60 MQStar accompanied by a 1 MQCompanion are shown 
in Fig. 2. A value for Oo smaller than reached at an earlier 
stage of the evolution implies that Roche lobe overflow 
would have occurred at that earlier moment. 

For each tabulated point of the evolutionary track we 
calculate the mass of the envelope Mg by subtracting the 
core mass from the total mass. Because the stellar evolu¬ 
tion models incorporate overshooting, which tends to in¬ 
crease the core mass, we calculate the core mass by mul¬ 
tiplying the value found from Eq. 2 with 1.125 (Maeder 
& Meynet 1989). We also know the minimum separation 
after spiral-in for a detached binary. From this we cal¬ 
culate the minimum separation at the onset of spiral-in 
with Eq. 4, and the minimum separation of the initial bi¬ 
nary ao,m with Eq. 6 . These minimum separations are also 
shown in Fig. 2. 

The value of Oo of the first evolutionary point at which 
Oo > ao,m, which we denote as amin, corresponds to the 
minimum initial separation that the binary must have to 
survive the spiral-in. The maximum of the values of Uq for 
all evolutionary stages, indicated with ao,max, corresponds 
to the maximum initial separation of the binary at which 
the primary can reach its Roche lobe. Only binaries with 
an initial separation in the range amin ~ o,o,max can evolve 
into low-mass X-ray binaries. For a 20 MQprimary with a 
1 M 0 secondary this range is 1000 - 1590 i?©. 

The Roche lobe can only be reached for the first time 
in those evolutionary stages for which Uo is larger than 
the Oo’s at all earlier evolutionary stages. Those stages are 
marked in Fig. 2 with •. Note that core hydrogen burning 
ends in tabulated point 13, and helium core burning begins 
in point 21. Helium core burning ends in point 43, and car¬ 
bon burning starts in point 46. For a star of ~ fiM©, the 
radius of the star expands following the end of core hydro¬ 
gen burning, and mass transfer during this first ascent of 
the giant branch is called case B. At the onset of core he¬ 
lium burning, the star shrinks. It expands once more after 
the end of core helium burning, and mass transfer during 
this second ascent of the giant branch is called case C. For 
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Fig. 2. A: • and -|- indicate radii for subsequent evolution¬ 
ary stages of the 20 Mq primary as tabulated in Schaller et al. 
(1992). From the mass of the primary at the tabulated point 
one may calculate the semi-major axis of a binary with a 1 Mq 
secondary in which the primary fills its Roche lobe, and from 
this the semi-major axis Uo, shown as a solid line, of the binary 
at the beginning of its evolution. For each tabulated evolution¬ 
ary stage we calculate the minimum semi-major axis at which 
the core of the primary survives the spiral-in of a 1 Mq com¬ 
panion, and from this the minimum semi-major axis Oo.m, indi¬ 
cated with the dashed line, of the binary at the beginning of its 
evolution. Primaries at the evolutionary stages marked with a 
-I- cannot fill their Roche lobe for the first time at that stage but 
will reach their Roche lobe at an earlier point in their evolution. 
A low mass X-ray binary is formed when Uo > ao,m for evolu¬ 
tionary stages indicated with •. B: as A for a 60 MQprimary 


the 20M© star shown in Fig. 2, however, the radius does 
not shrink at the onset of each new phase of core fusion, 
but continues its expansion throughout its evolution, once 
the Hertzsprung gap is passed. The 60M© shrinks at the 
onset of helium fusion in the core, mainly due to extensive 
mass loss. 

As shown by Fig. 2 a 60 M© primary with a 1 M© sec¬ 
ondary can, according to the same reasoning, only evolve 
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into a low-mass X-ray binary if its semi-major axis is in 
the very small range of 980 - 1100 Rq. 

We determine the mass and size of primaries in a 
range of masses at the moment they fill their Roche lobes 
at ao,max for each tabulated stellar evolution track. The 
masses and radii at ao,max of the stars that are not tabu¬ 
lated by Schaller et al. (1992) are determined by a linear 
interpolation between the tabulated models. The resulting 
values for ao,max are shown as a solid line in Fig. 3. 

If Roche-lobe overflow for all initial semi-major axes 
smaller than ao,max leads to a merger, then amin is not 
properly defined. A lower limit can be obtained by com¬ 
puting amin with the stellar parameters that correspond 
to the point where ao,max is reached. 

Mass loss from stars more massive than ~ IOOMq is so 
copious that these stars lose their entire hydrogen envelope 
before they expand on the giant branch, i.e. before Roche- 
lobe contact is achieved. The common-envelope phase 
hardly leads to a spiral-in, whereas further attrition of the 
core to a small final mass (~ 8 Mq according to Schaller 
et al. 1992) causes the binary orbit to expand. As a result 
the final orbit is so wide that the 1 Mq secondary never 
reaches its Roche lobe. 
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Fig. 3. Lower limit to the initial semi-major axis at which the 
binary survives the spiral-in, as function of the initial mass of 
the primary, calculated with use of the evolutionary sequences 
by Schaller et al. (1992). The upper (lower) dashed line gives 
the limit determined from the condition that the secondary 
star (helium core) is smaller than its Roche-lobe, The solid line 
gives the upper limit at which the primary reaches its Roche 
lobe at its maximum radius of Rmax = lOOOR©. The secondary 
is assumed to be a 1 M 0 Star. The dotted line indicates the 
initial semi-major axis of a binary in which a 1 Mq secondary 
can fill its Roche lobe after the primary has lost its entire 
envelope without ever having filled its Roche lobe 


Figure 3 shows, as a function of the zero-age mass of 
the primary, the lower limits to the semi-major axis of the 
initial binary at which the binary survives the spiral-in 


Umin (or its Conservative lower limit), and the upper limit 
for which the primary reaches its Roche lobe ao,max- 

Comparison of Fig. 3 with Fig. 1 illustrates that the 
mass loss of massive stars and the concurrent widening 
of the binary reduces the maximum initial separation for 
which the binary reaches Roche-lobe contact. For stars 
with ~ 70 Mq the maximum stellar radius is smaller than 
1000 Rq, which leads to a further reduction of ao,max- The 
wind mass-loss affects the minimum initial separation nec¬ 
essary to survive the spiral-in in two ways: the widening 
of the orbit reduces this separation, whereas the reduction 
of the envelope mass enlarges it. 

Figure 3 indicates that the formation rate of low-mass 
X-ray binaries with a neutron star greatly exceeds the 
formation-rate of low-mass X-ray binaries with a black 
hole, because the range of allowed initial separations is 
larger for neutron star progenitors. 

Romani (1992) mentions the possibility that the sec¬ 
ondary star might fill its Roche lobe as it evolves on the 
asymptotic giant-branch after the primary collapsed into a 
remnant without ever having reached Roche-lobe contact. 
If this happens the binary orbit also widens dramatically, 
according to Eq. 5. If the orbit widens too much, the sec¬ 
ondary will not reach its Roche lobe. In Fig. 3 we show the 
maximum initial semi-major axis for which a 1 Mq star 
can reach its Roche lobe after its companion has lost its 
envelope. It is seen that this semi-major axis is so small, 
that it invalidates the assumption that Roche lobe contact 
of the primary has been avoided. 

4. Galactic formation rates 

The formation rate of binaries that reach Roche-lobe con¬ 
tact and survive the spiral-in can be estimated by com¬ 
puting the birth rate of binaries with suitable initial pa¬ 
rameters; primary mass, mass ratio and semi-major axis. 
For simplicity we neglect the eccentricity. For the initial 
mass function for the primary we use a Salpeter function 
(Salpeter 1964) integrated over the Galaxy: 

Tg(M) = 0.05M-2-35 [M© yr-i]. (7) 

This normalization gives a supernova rate in agreement 
with the observed rate. The initial semi-major axis distri¬ 
bution r(a) is taken flat in logo (Kraicheva et al. 1978). 
We assume a flat initial mass-ratio distribution, $(( 7 ) = 1. 
For a given mass of the primary M the fraction of binaries 
with a mass of the secondary between (m — e) and (m -I- e) 
is then proportional to 2 e/M. 

Primaries with an initial mass between 10 and 40 M© 
are assumed to leave a neutron star after the supernova, 
progenitors with a mass between 40 and 100 M© form 
black holes. The mass of the secondary is taken to be 
1 Mq. The minima and maxima for the initial semi-major 
axis are computed as described in the previous section. We 
integrate the initial distribution functions for the primary 
mass, the mass ratio and the semi-major axis between the 
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mass limits Mmin and Mmax and the corresponding limits 
for the semi-major axis Omm and ao,max (see Fig. 3): 

P Mmax pO-o,max 

Br= / 'i/G{M)—r{a)dadM[yr~^]. (8) 

JMmi„ Jamin 

With e = 0.15 Mq, with a binary fraction of 50% of high 
mass stars, and with the lower limit to amin set by the 
main-sequence star, the galactic formation rate of binaries 
that reach Roche-lobe contact and survive the spiral-in is 
2.2 • 10“®yr“^ for the binaries with a neutron star and 
9.6 • 10“®yr“^ for the black-hole binaries. The formation 
rate thus calculated for the low-mass X-ray binaries with a 
neutron star is compatible (within the rather wide uncer¬ 
tainties) with the rate derived from the observed numbers 
of such binaries. The formation rate of low-mass X-ray bi¬ 
naries with a black hole is about 1% of the formation rate 
of low-mass X-ray binaries with a neutron star, whereas 
the observations indicate equal formation rates for these 
two types of binaries. 

5. Discussion 

We briefly discuss the assumptions and model uncertain¬ 
ties that affect our computations most. 

In our calculations the effect of convective overshoot¬ 
ing on the mass of the helium core is taken into account 
in a relatively simple manner, since the tabulated stellar 
evolution tracks do not provide masses for the stellar he¬ 
lium core. The importance of convective overshooting is 
uncertain, and therefore we investigate small variations in 
its effects. 

If we increase the mass of the core by increasing the 
multiplication factor for Eq. 2 discussed in paragraph 2 of 
Sect. 3 from 1.125 to 1.250, the lower limit amin to the 
semi-major axis for which a binary survives the spiral-in 
is slightly reduced, and consequently the birthrate of neu¬ 
tron star binaries is enhanced by ~ 28%, and the birthrate 
of black hole binaries increases with a factor ^ 4. An in¬ 
crease of the exponent in Eq. 2 to 1.57, which is not un¬ 
realistic according to Iben & Tutukov (1985), reduces the 
discrepancy between the birthrates for black holes and 
neutron stars in low-mass X-ray binaries to a factor 15. 
Thus neither of these changes in Eq. 2 leads to a suf¬ 
ficiently high birthrate of black-hole binaries relative to 
neutron-star binaries. 

The efficiency a\ at which the common-envelope is 
expelled upon the spiral-in is highly uncertain. The effect 
of varying aA on the galactic formation rate for binaries 
that reach Roche-lobe contact and survive the spiral-in 
is demonstrated in Fig. 4. Decreasing aA results, as ex¬ 
pected, in a strong reduction in the formation rate of 
low-mass X-ray binaries with a neutron star as well as 
with a black hole. By increasing aA to a value larger than 
about unity, the formation rate of low-mass X-ray bina¬ 
ries approaches a constant rate, limited by the number of 


primordial binaries formed in the galaxy. Note that aA 
might have different values for neutron star- and black- 
hole progenitors. 

To investigate the effect of the stellar metallicity on the 
formation rate of low-mass X-ray binaries we compute two 
more models using the evolutionary tracks calculated by 
Charbonnel et al. (1993) for Z = 0.004 and Schaerer et al. 
(1993) for Z = 0.04 (see Fig. 4). Apart from the enhanced 
mass-loss rates in the model with high metallicity which 
results in a larger formation rate of low-mass X-ray bi¬ 
naries, these model variations have little influence on our 
conclusions. The models with a lower metallicity lose less 
mass in the stellar wind and the fraction of binaries that 
survive the spiral-in decreases accordingly. 



Fig. 4. Birthrate (per year in the galaxy) of low-mass X-ray 
binaries with a neutron star (upper lines) and a black hole 
companion (lower lines) as function of the common-envelope 
efficiency-parameter a A. The full lines are from the stellar evo¬ 
lution models from Schaller et al. (1992), the dotted lines are 
from the low metallicity models with Z = 0.004 from Char¬ 
bonnel et al. (1993) and the dashed lines are computed with 
Z = 0.040 from Schaerer et al. (1993) 

The limit of ~ 4OM0, which we use throughout the 
calculations, was originally derived from the observation 
of the high-mass X-ray binary LMC X-3 (van den Heuvel 
& Habets 1984). Recently the analysis of the data of 
Wray 977, the Bla hypergiant companion of the X-ray 
pulsar GX301-2, has increased the empirical lower mass 
limit for the formation of a black-hole in a binary to 
Mo ~ 50 Mq (Kaper et al. 1995). This only exacerbates 
the problem of too low a formation rate of black-hole bina¬ 
ries in our calculations. Maeder (1992) argues on the basis 
of the galactic ratio of the abundances of metals and of he¬ 
lium that the lower limit for the formation of a black hole 
is 20 - 25 Mq. As can be seen in Fig. 3, such decrease of 
the lower mass-limit has little effect on our computations, 
because amin > ao^max in the range of 25-40 Mq. 

The stellar evolutionary tracks of Schaller et al. (1992) 
follow the single stars beyond the point of carbon burning. 
The large amount of mass and angular-momentum lost 






6 


S.F. Portegies Zwart, F. Verbunt & E. Ergma 


in the stellar wind during earlier phases of the evolution 
of the more massive primaries cause the binary orbit to 
widen and prohibits mass transfer afte the onset of core 
helium burning (see Fig. 2). For the amount of angular 
momentum lost per unit mass in the stellar wind we used 
the standard description of isotropic mass loss from the 
donor (see Eq. 5). Mass transfer at an evolutionary stage 
well after the onset of helium fusion (say, beyond point 25 
in Fig. 2) can only be an effective channel for the forma¬ 
tion of low-mass X-ray binaries with a neutron star as the 
compact object, if the loss of angular momentum is high 
enough for the orbit to shrink. For a more massive pro¬ 
genitor, larger than ~ dOM©, which leaves a black hole 
instead of a neutron star after the supernova, such late 
mass transfer is not likely to become important even for 
such high loss of angular momentum: the star has evolved 
into a naked helium core and a merger will be unavoidable 
upon Roche-lobe contact. 

The evolution of the helium star and the effect of the 
supernova on the binary orbit are neglected in our calcula¬ 
tions: both effects tend to lower the formation rate of low- 
mass X-ray binaries (see e.g. Portegies Zwart & Verbunt 
1996). Mass transfer from the helium star to its companion 
can result in a merger and the asymmetry of a supernova 
can dissociate the binary. Both effects, however, tend to 
reduce the formation-rate of neutron stars as well as black 
holes in low-mass X-ray binaries and do consequently not 
solve the birth-rate problem discussed in this paper. 

6. Conclusions 

Our computations, in contradiction to the results of Ro¬ 
mani (1992), predict a formation rate for low mass X-ray 
binaries with a black hole which is much smaller than the 
value derived from the observed numbers and estimated 
X-ray lifetime. The striking difference between our and his 
results can be seen immediately by comparing his Fig. 1 
with our Fig. 1. Romani’s values for Umax are much higher 
than ours: this is due to the fact that he uses l/q instead 
of q in the equation for the Roche lobe (our Eq. 1). (This 
error has been silently corrected in Romani 1994.) This 
is true both for binaries in which mass transfer preceeds 
the supernova explosion, and for binaries in which the pri¬ 
mary explodes after losing its envelope without ever hav¬ 
ing filled its Roche-lobe. The discrepancy between theo¬ 
retical and observed formation rates cannot be solved by 
invoking different metallicities for the progenitor systems, 
nor by assuming different efficiencies for the envelope ejec¬ 
tion during spiral-in. 

Black-hole binaries can be produced in larger numbers 
only if it is assumed that stars with initial masses less than 
approximately 20 Mq can collapse to black holes; or alter¬ 
natively if it is assumed that the angular momentum loss 
caused by the stellar wind is so high that the binary orbit 
shrinks; or alternatively if the collapse of a helium core to 
a black hole is asymmetric, so that the post-supernova or¬ 


bit can be smaller than the pre-supernova orbit (see Porte¬ 
gies Zwart & Verbunt 1996). 
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